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Abstract 

Three-point vertex diagram plays a key role in the whole renormaliza- 
tion program of several QFT (quantum field theory) models such as QED, 
QCD, the Standard Model of eletroweak interactions and so forth. The 
exact analytic result for the triangle diagram therefore is fundamental. 

In this work we calculate in two different ways a two-point two-loop 
massless Feynman diagram using what we call a "gluing" technique in 
the context of NDIM (Negative Dimensional Integration Method). The 
two-loop diagram in question can be "glued" in two different ways and 
we show that both yield the same result and reproduce the one calculated 
via NDIM for the complete diagram, which, of course, is equivalent to the 
exact solution obtained by normal positive dimensional calculation. 

Furthermore, in the process we conclude that the usual massless off- 
shell triangle diagram result does not hold anymore and present a new 
solution for it with only three hypergeometric functions F4. 
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1 Introduction 

Today in quantum field theory, it is necessary to calculate increasingly complex 
Feynman diagrams as the theory and the experiments require a higher accu- 
racy of the scattering amplitudes; a very good example of this being Kinoshita's 
quest: to calculate (g — 2) up to a 5 order pQ. Several techniques have been 
applied for that purpose — most of them in the context of dimensional regu- 
larization [5] or analytic regularization [5] — and among them we can mention 
the powerful Mellin-Barnes contour integration [4j [5J, [6] , the method of Gegen- 
bauer polynomials [7J, the differential equations technique [5] and others [S]. 
The NDIM developed by Halliday and Ricotta [TO] has shown itself as a reli- 
able one when applied to the calculation of diagrams of one- [ITJ [T3] , two- 
and multi-loops [141 115] , with scalar and tensorial structures and in noncovari- 
ant gauges [16] . One of the advantages of NDIM is that it allows us to avoid 
the often cumbersome parametric integrals, transfering the problem into easier 
solving systems of linear equations instead. Another advantage of NDIM is that 
the exponents of propagators are taken to be arbitrary integers, so that one can 
solve the general case for each type of graph. 

Despite of having these advantageous features which turn itself an efficient 
and simple method, NDIM has some drawbacks pointed out in earlier works: the 
large number of systems of equations and the difficulty in dealing with the ever 
increasing complexity of the hypergeometric type series that results. Related to 
this difficulty is, for example, the summing up analytically of hypergeometric 
series with unit argument, such as, p F q (...\l) that often appears in two-point 
function calculations. The approach to overcome the first difficulty, namely, to 
reduce the growing number of linear systems — it grows with the number of 
loops and legs attached to the diagram at hand — was presented by Gonzalez 
and Schmidt for massless diagrams: they proposed a way to write down the 
generating gaussian integral in order to optimize (to minimize!) such number. 
Then the whole calculation can be made simpler and faster, and the number of 
hypergeometric functions p F q (...\l) left in the final result is also minimum. 

A great effort is being conducted in order to study maximally supersym- 
metric Yang-Mills theory (MSYM). Several authors [17] are tackling a rather 
difficult task: to calculate certain scattering amplitudes exactly since Maldacena 
conjectures that higher loop contributions can be written in terms of one-loop 
amplitudes. Tests of this conjecture have been realized from 4-point 2-loops 
to even more challenging 4-point 5-loops Feynman amplitudes, involving the 
so-called dual conformal integrals. The outcome of this program can be a re- 
summation of the entire perturbative series for a given physical process. The 
reason to have powerful methods to tame these integrals is clear. 

We present in this paper another way to apply NDIM in Feynman integral 
calculations, the "gluing" approach. In general a Feynman diagram, e.g., the 
2-loop master diagram, is represented by an integral, 
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where p is the external momentum, and i, j, l,m,n are exponents of propagators, 
which can be made arbitrary in the whole calculation. One could rewrite the 
above integral as, 



and we readily recognize the integral in momentum q as an off-shell triangle 
one, which has a well-known result 18J that can be written in terms of four 
Appel hypergeometric functions of two variables F±{. . .\k 2 j 'p 2 , (k — p) 2 /p 2 )- It 
is straightforward to see that the remaining integral is a self-energy one with 
shifted exponents of propagators, 



where T is a factor which depends on p 2 as well as the exponents of propagators 
and dimension D. The new exponents v\ and V2 also depend of the former ones 
i,j,l,m,n, as well as of dimension and sum indices of F4 functions. However, 
straightforward application of this does not yield the correct result. Here comes 
an important point: to carry out the second integral ([3]) one has to perform the 
integral over the whole space, for this reason the result of the former one must 
hold on the whole range of momentum k. The well-known result of the off-shell 
triangle, written as a sum of four Appel's hypergeometric functions Fi{...\x,y), 
is not valid for every momentum; these momenta must be such that \x\ < 1, 
\y\ < 1 and \y/x\ + \y/y\ < 1. In other words, the series is defined inside some 
region of convergence and for this reason the well-known result of Boos and 
Davydychev [T3J, [TS] can not be used in fl2J) . 

In this paper we use NDIM to solve a massless two-loop self energy diagram 
in a different approach as used before [19j . namely, integrating loop by loop. 
We separate the diagram into two simpler parts, each one a single loop diagram 
itself and then solve these parts "gluing" them to obtain the final result. This 
technique of dealing with subdiagrams could simplify the solution of larger dia- 
grams that would lead to difficult systems with a great number of equations and 
variables. As far as we know only Bierembaum and Weinzierl [20] studied some 
Feynman diagrams using similar ideas and Mellin-Barnes method. Kostrykin 
and Schrader presented a generalized Kirchoff rule to "gluing" quantum graphs 
and they relied on a complicated star product [21]. "Gluing" diagrams is not 
a trivial task in quantum mechanics nor in quantum field theory. However, we 
think that if implemented such method can simplify computations, being very 
easy and reliable, since in order to calculate a multi-loop diagram one could, 
in principle, use only well-known one-loop integrals. The main objective of the 
present work is to elaborate this program within the NDIM context. 

There are two manners to cut the two-loop diagram — here identified as 
"flying saucer" diagram — in the particular case presented in this work, as 
will be shown in Sections 3 and 4. The first manner involves a one-loop self- 
energy diagram that we call "sword-fish" , which is very simple to calculate. The 
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second one uses the one-loop off-shell triangle diagram, which is written in terms 
of only three hypergeometric functions F4 — the form suitable for integrating 
the second loop — , a simplified result accomplished by invoking momentum 
conservation [55]. 

The outline of our paper is as follows: in Section 2 we show the integral 
pertaining to the "flying saucer" diagram and present the approaches to make 
the cutting. In Section 3 we solve first the "sword-fish" diagram, whereas in 
Section 4 we do it via the one-loop triangle diagram first and in Section 5 we 
present our concluding remarks. 



2 The flying saucer diagram 

We will consider a particular case of the flying saucer diagram, the so called 
"side view", as shown in Fig. 1. There is another type of this graph, which we 
call flying saucer "front view" . The difference between them is the exponent of 
the propagator of the particle with momentum k. Clearly the following gaussian 
integral is related to this diagram, 

I(p 2 ;D) = J d D kd D qexp[-aq 2 - [3k 2 - 7 (p - k) 2 - u{k - q)\ (4) 

in other words it generates the Feynman loop integrals, see eq.©. 
One can integrate it and compare it with its own Taylor expansion, 

I(p 2 ;D) = £ (-l) i+j+l+m ^^^id,l,rn-,D-,p% (5) 

i,j,l,m=0 



where 

-s( l ,j,l,m;D;p 2 )= I d D kd D q(q 2 y(k 2 y[(p - k) 2 ] l [(k - q) 2 ] m , (6) 



which is our complete negative dimensional integral. We will show all the nec- 
essary steps for the cut cases in more details in the next sections, but for now 
it is enough to know that the solution of the complete diagram is, before the 
analytic continuation [19j . 

aftAi.m^-M'W r ( i + i >r(i + 0r<i + m) 



r(l - a - \D)T(l + i + j+m + |D)r(l —i — m—D 



r(i + <r)r(i - i - |D)r(i - m - \d) 

(7) 



r(l - l - \D)T{1 + i + m + |-D)r(l + l-a) 



where a = i + j + I + m + D. 

At this point, our new proposal is to cut the diagram, solve each part sep- 
arately and then "glue" them together. We expect, of course, that the final 
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p-k p-k 

Figure 1: Two-point two- loop Feynman diagrams: Fying-saucer side-view and 
front-view. 

result should be the same as in ||7J). In Fig. 2, we show the two manners to cut 
this diagram. The first way leaves a one-loop sword-fish graph plus a vertex, 
whereas the second way leaves a one-loop triangle and a propagator. In the next 
sections, we solve them in details. 

3 The sword-fish diagram 

We begin with the generating integral ((4|), where one can integrate first in the 
q momentum, 

I(p 2 ;D)= J d D kexp[-pk 2 -7(p- kf] J d D qexp[~aq 2 - uj{k - q) 2 }, (8) 

which represents a one-loop with two massless propagators. Completing the 
square, the q integral can be solved easily, giving 



I q (k i ;D)=[——) exp 
a + lo ' 



otujk 



a + u) 



(9) 



Expanding this result in Taylor series and using the multinomial expansion, 
one obtains 

I q {k 2 - D) = ( ] ( m 2 ,^ , (10) 

^— ' ni!n2!n 3 ! 

with the constraint — ri\— — 712+713 coming from the multinomial expansion. 
The next step is to also expand the q gaussian integral ([5]) in Taylor series, 

oo i m 

I q (k 2 ;D)= £(-l)W^G g (i,m;Afc 2 ), (11) 

h3=0 

where 

m; D; fc 2 ) = / ^(^[(fc - g ) 2 ]"\ (12) 
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k - ci 



p-k 



Figure 2: Two-point two- loop Feynman diagrams solved using a "gluing" tech- 
nique. We show two ways of cutting them in order to integrate loop-by-loop. 



This is the main sword-fish integral in negative dimension. The exponents 
of the propagators are positive and one has to consider a negative D. In the end 
of the calculation, one does an analytic continuation of the result to negative 
exponents and positive D. Comparing (|10p and (|lip . an expression for the 
negative- 1? integral can be written, 

Q q (i,m;D;k 2 ) = (— k)^ V] — i — i — Ani+n 2 <Wi+n 3 > ( 13 ) 

f— 'ni!n 2 !n3! 

rti— 

where a product of gamma functions is defined, 

G = T(l + i)T(l + m)r(l + i + m - 2m). (14) 

The Kronecker's deltas in (j!3p and the constraint of the multinomial expan- 
sion of pop form a system of three equations and three unknowns, 

ni + ri2 = i 

n% + ria = 77i (15) 

77l + 71 2 + 77-3 = - J-D 

that can be easily solved with a unique solution. The integral ^s q (i, m; -D; fc 2 ) is 
then given by 

9f,(t, m; D; k 2 ) = (-n)^ (k 2 ) t+m+ ^ D P(i, m; D), (16) 

where 

r(i + i)r(i + m)r(i-i-m-i>) 

P(i,m;D) = — : — - — : — - — i-^r- (!') 

v ; r(l + i + m + i£>)r(l-i-iD)r(l-m-i£») 

With this result, it is straightforward to find the solution of the complete 
flying saucer diagram. From the Taylor series ([6|) , it can be seen that the integral 
in q is already done, so putting (fTS]) in ©, one has 

^j,l,m;D;p 2 ) = P(i,m;D) f d D k(k 2 y+*+ m +^ D [(p ~ k) 2 } 1 , (18) 
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but this is exactly the integral (fl"2")l with the variables *Sk(i+j+m+-hD, I; D;p 2 ), 
that is, 

j, l,m;D;p 2 ) = (-tt) d (p 2 y+^ l+m+D P(t, m; D) 

x P(i+j + m+~D,l;D), (19) 

and, finally, 

%(i,j,l,m;D;p' 1 ) = (-Ti)(p z ) 
r(l - a - \D)Y{1 + i + j + m + ^D)T{1 -i-m-D 



2 , , 2 , ff r(i + i)r(i + or(i + m) 



r(i + cr)r(l - i - 2-D)r(i - m - ±D) 

(20) 



T[l-l-\D)T(l + i + m + \D)T{l + l-(x) 

which is exactly ([7]), with a — i+j + l + m + D. This procedure to solve 
the flying saucer diagram is much simpler and faster than solving the complete 
diagram as it was done in 19J. 



4 The off-shell one-loop triangle diagram 

In this section we present the other mode of cutting the complete diagram. 
This is by far the most difficult and laborious way; however considering the 
importance of the triangle diagram and the new result that we will show justify 
the whole process. Instead of ([8]), one could begin with the k integral: 

I(p 2 ;D) = J d D qexp[-aq 2 } J d D fcexp[-/3fc 2 - j(p - k) 2 - u;(k - qf], (21) 

that corresponds to a one-loop three point function. Working on the k integral, 
it gives without trouble, 



I k (p 2 ,q 2 ;D) = 



exp 



(3jp 2 + /3ojq 2 + yur 

(3 + J + LU 



2 1 



with r = q — p. Expanding (|22p in Taylor series, one gets 

(-l)"i+ n =+™3(_ ni - n 2 - 713 - \D)\ 



I k (p 2 ,q 2 -D) = ^)- 



x (p 2 ) 711 (q 2 ) 112 (r 2 )™ 3 / g™ 1 +™ 2 +™ 4 ^ ni + n3 +" 5 [J j™ 2 +™ 3 +™6 



(22) 



(23) 



with the constraint ri\ + 112 + Tis + + n.5 + n@ = —\D coming from the 
multinomial expansion. Now one expands the original k integral (|2ip in Taylor 
series, obtaining 



2 ,,:>. /-)• — V tl ft]^, 1 ,.,™^, I/, I /I..,-' .,2, 

j\l\m\ 



I k (p 2 ,q 2 ;D)= Y, 

j,l,m=Q 



-^ l cj m Q k (j,l,m;D;p 2 ,q 2 ), (24) 
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where the corresponding integral in NDIM is 



Comparing 
general relation 



l,m;D;p 2 ,q 2 ) = I d D k(k 2 y[(p - k) 2 ] l [(k - q) 2 ] m . (25) 
and (|2"l)l by its /3, 7 and ui powers, the integral 3^ has a 



9f fc (i, l,m;D-p 2 iq 2 ) = (^(-l)"^"™ 

x b 2 )" 1 ( q 2 r (r 2 r5 3 , ll+n2+n ^ ni+v S ,„.S 



E 



(-1) 



G 



ni!n2!n3!n4!n5!ng! 



where 



G = r(l + j)T(l + l)T(l + m)r(l - m - n 2 - n 3 - -£>). 



(26) 



(27) 



Considering the deltas in (|2"6"|) and the multinomial expansion in (|2"3")l . all the 

constraints of the problem are 



Hi + n 2 + n 4 = j 
n\ + n 3 + rt 5 = I 
n 2 + n 3 + n 6 = m 
n\ +n 2 +n 3 = 4> 



(28) 



where <f> = j + I + m + \D. 

Therefore one has a system of four equations and six variables. There are 
C\ = 15 possible ways to solve the system leaving two free variables, that 
is, ending up with a double series. From these 15 possibilities, 3 have zero 
determinant, so there are only 12 non-trivial solutions that can be grouped 
together according to the ratio of the momenta, or in other words, according to 
the kinematical configuration. Usually, the groups are made by four solutions 
and each one is linked to the others considering the symmetries of the diagram. 
As said before, the leftover variables of the sum form a double series that can 
be written in terms of Appel hypergeometric functions F4 [2 3) , 



oc 



Fi (a, b, c, d; x, y) 



(a) m +n(b) m+n x m y 

m,n—0 



(c)m(d) 

with the Pochhammer symbol designated by 

T(a + to) 



m\n 



(a\m) = (a) r 



r(a) 



(29) 



(30) 



and obeying the useful relations 

(a\m + n) = (a\m,)(a + m\n) 



(a\ 



TO 



LIT 

(1 — a j to) 



(31) 
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In terms of these definitions, the first set of solutions is 



l,m;D;p 2 ,q 2 ) = (-ir)- 
x A X F 4 ( -j, -<j>,l + 1 - <j>,l + m ■ 



2 2 

q p 



+A 2 F 4 [-1, -j - I + </>, 1 - I + 4>, 1 + m 



2 2 
9 P 



+A3F4 —m, —j — m + 4>, 1 + / — 4>, 1— m + <f>: 



2 2 
<7 P 



\-l-m 



1 



D, 1 - 1 + 6, 1 - m- 



<? 2 P 2 



(32) 



where the multiplicative factors are 
At 



2^ (l + 0|-20-^D) 
r J (l+Z|-0)(l+m|-0) 



-4o 



-4, 



(* 2 )*(S)" 
(P 2 )'(^) 



(i+0-;|-20+i-i_D) 

(l+j|2-0)(l+m|-0) 

1 (l+0-m|-20+m-il)) 
(l+j|m-0)(l+Z|-0) 



(33) 



(*H*) 



-/ 



(l-j-lD|2j+lg) 
(l+/|i+iZ5)(l+m|j+i_D) 



This result agrees with the one calculated in the first reference in [4]. The 
other two sets with the remaining eight solutions can be found making the 
replacements 



%(j,hrn;D;p 2 ,q 2 ) = %(q^r,j 
$l(j,l,m;D;p 2 ,q 2 ) = Qf£(p <-> r,j 





m) 



(34) 



Now, it comes the crucial new step that will permit us to complete the inte- 
gral (|2Tj) . As presented in (|32[) . the solution is not valid in the whole space, nor 
takes into account the momentum conservation r = q — p. This last constraint 
subtly reduces the number of functions F4 from four to three because those orig- 
inal four are not linearly independent as they should be to form a basis. Two 
of them can be rewritten in terms of another function F4 that pertains to one 
of the other sets of different kinematical regions (f34|) . So, with the constraint 
of momentum conservation, there are only three linearly independent hyperge- 
ometric solutions F4 and they hold on the whole space. In order to reduce the 
four functions F4 in the solution (132j) , we need to combine two of them to give 
another F4 in a different kinematical region. The right way to do this combining 
is given in [22] considering the analogy between Feynman diagrams and electric 
circuits. Then, applying this to our case, we need to keep the first and the third 
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terms, and the second and fourth ones are combined together using the relation 



F 4 (a, b, c, d;x,y) 
( a )-fc 



(<*)-. 



X 1 



(-y)~ a .F4 a, a + 1 - <2, c, a + 1 - 6; -, - + 



y y 



+ hr^(~y) F * [b+l-d,b,c,b + l-a; 



x 1 



(35) 



(d)-b \ y y, 

to see that the resulting F4 is exactly a function that appears in 3^(j, I, m; D;p 2 , q 2 ) 
given in (f3"l)) . So (|32p is in fact given by 

^(j,Z,m;£>;p 2 ,cZ 2 ) = (-tt)*x 



2 2 

. 1 j F: j -j, -<f>, 1 + J - 0, 1 + m - & ^L, ^ 



( -m, -j - m + 0, 1 + / - 0, 1 - m + 0; — , — 



2 2 

-A5F4 ( -l,-m-l + 4>,l+j-<f>,l-l + <f>\^ s ,^ s 



(36) 



where 



A 5 = (q 2 Y 



p 2 ) (l + m\l-<t>)(l+j\ -< 



The kinematic region of the last term is distinct from the previous ones, 
but as can be seen in Fig. 3, these two regions are connected, allowing the 
momentum to hold from — oo to +00. 

With the solution (|3"6")) in hands, it is now possible to finish the work. Going 
back to ©, the k integral is done, so putting ([3"6) in it, 



S(i,j,l,m;D;p*) = I d D q(q 2 y%l(j,l,m;D;p',q'), 



(37) 



which gives us three integrals of self-energy type, e.g. |T! 

^,3,l,m;D;p 2 ) = (~n)% £ / d D ql A[(q 2 y+* \(q - pf 



,6=0 



-a—b 



m — a — b 



+A' 3 (q 2 Y+ a {q-pf + A 5 {q 2 f +l - l+a (q - pf 



(38) 



where 



At 
A 3 
A' 



HUNW (1 + g-ag-jg (p 2 ) b 
(1 + j - 0) a (l + m- <t>) b (1 + 0-0 (1 + ™)_ 4> a\b\ 
(-m) a+b (-j - m + 0) a+b (1 + - m)_ 20+m _ i D (p3) ^+^ m 
(1 + Z - 0) o (l - m + 0) b (1 + j) m _, (1 + I), a\b\ 



(-l)a+b(-m - I + (f)a+b (1 + <f> - 0- 2 <A-H-ii3 (p 2 ) 
(1 + j - 0) a (l - j + 0) b (1 + m)^ (1 + a\b\ 



I — a — b 



(39) 
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Figure 3: Region of convergence of the hypergeometric series F4 when one does 
the correct transformation in order to account the momentum constraint. 



Using the solution (|16p with the right replacements of the variables, one gets 



S(i,jJ,m;D;p 2 ) = h + I 3 + 1 5 , 



where 



>-2<j>-iD 



(l + O_ (l + m)_ 

r(i + z)r(i + 0)r(i -i-4>-p) 
" r(i + i + <t> + |D)r(i - i - |i>)r(i -<i>-\d) 

^ (~i)»+6(l + *)a(* + |£>)a(l - t - - D) 6 (-» 

: L — 

a,b=0 



\D) b 



o!W(l + J - 0) a (l + "1 - 0) b (l - - ±£) a+b 



(40) 



(41) 



h = (-*•)* (P 2 ) 



(-iy 



(l + ^)-2^-iD(- 



,(1 + 1) + 



T(l + i)T(l + m)T(l-i-m-D) 
r(i + i + m + i£))r(l - i - ±£>)r(i - m - |D) a ^- o 



(l + *)a 



a!6! 



(1 - i - m - D) 6 (i + \D) a {-i TTL |D) 6 (-j - TO + 4>)a+b 



(1+ I - 0)a(l - "1 + 0) h (l - m - |D)a+6 



(42) 
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h = (-7r)~ [P 

r(i)r(i 



2\ i+4>+iD 



(1 



-20- 



(- 



(l+i)-*(l + "O-*(l + "»-0i 
{l) b (l + i-l + cP)t 



' T(l + i - I + <j> + |-D)r(l - \D) 



E 

a,b=0 



a\b\{l + i-l + ^+\D) a+b 



(i-l + <j>+ \D) a {\D) b {-l) a+b {-l - m 



)a+b 



(i -i- +)„(! + j - 0) 6 (i - 1 + 4> + D) a 



+b 



(43) 



The summations above can not be written in terms of a known hypergeo- 
metric series with defined properties. So one has to specify the values of the 
indices i,j,l,m and make the analytic continuation. Using (|31[) . it can be seen 
that due to the analytic continuation of the terms 



r(i) 



in (|43j) . I 5 vanishes. 



Making i = — 1 in I\ and 23, the terms (1 
index a and one has 



4 \ _ EM 
'/« — r(o) 



r(i-i£>) 

cancel the sum in the 



h = (-7T) (J3 ) 



(1 



-20-i_D 



(l + l)_^(l+m)_, 

r(i + i)r(i + <^)r(i — i — <f> — D) 



v{i + i + 4 

X 2^1 (-i, 1 



i£>)r(l - i - \D)T{1 -4>-\D) 



DA 



2 
m 



0;i) 



(44) 



(i + j')- (i + O- (i+i-^) m 

T(l + i)r(l + m)r(l -i-m-D) 
x r(i + i + m + fojr(i - i - i£>)r(l - m - i£>) 
x 2-F1 (1 — « — m — £>, -j - m + 0, 1 — m + (j>; 1) , (45) 

where 

(a )m(fr)m z r ' 
(c) m m! 



T? 1 1, \ sr^ { a )m(b)m x , 

2 F X {a, b, c;x)= — n r ( 46 ) 



m— 



is a hypergeometric series that can be summed up when it has unit argument, 

1 (c — a) 1 (c — 0) 

Using the property (|4T[) in (|4~4l |4"5")) , making the analytic continuation and 
specifying i = j = l = m = — 1, the integral (|4D)) becomes 

3 AC ( l ,j,Z,m; J D;p 2 )-( 7 r) I3 (p 2 ) 1J 4 Va - ' 1 -i X 



r (3 - 1£>) r (D - 2) 

r (2 - \d) r (4 - l>) 



r(|£»-4) 



(48) 
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which gives the exact result of the flying saucer diagram if one makes the analytic 
continuation and the specification of the exponents in the solution ([7]). 

5 Conclusion 

We presented in this work a technique to integrate multi-loop Feynman integrals. 
In our approach one can carry out the integrals loop-by-loop using well-known 
results of straightforward one-loop diagrams. We point out the advantage that 
each diagram when properly calculated — with off-shell external legs and re- 
specting the constraint on the external momenta — can be used to solve an even 
harder diagram and so forth. In this approach each diagram can be considered 
as a building block of another one, with more legs and/or loops. The same 
method can also be applied to the two-loop master integral and a result valid 
for arbitrary exponents of propagators can be obtained. 
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